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1. INTRODUCTION 


A leading tool in the manufacturing process is the control chart 
that is applied to watch manufacturing shifts. A slight change in the 
quality target may cause big losses for a company. The variable con- 
trol charts are designed to monitor the process using measurement 
data. Nonconforming items can be monitored using an attribute 
control chart. According to [1], “the traditional Shewhart np con- 
trol charts are the statistical control scheme most commonly used 
for monitoring the number of non-conforming items.” For more 
details, the reader may refer to [1-4]. 


Usually, control charts work when data is derived from a nor- 
mal process. When data moves away from normality, the chart 
designed with normal distribution cannot monitor the process. 
Control charts designed with non-normal distributions present a 
good alternative. According to [5], “In most life testing, to reduce 
the test time of the experiment, a failure-censored (type-II) scheme, 
or time-censored (type-I) scheme is usually adopted” Therefore, 
the control chart designed with a time-truncated control can be 
used to save monitoring time. In the operational procedure of time- 
truncated charts, the experiment and watching times are fixed in 
advance, where an item is labeled as defective if its failure/service 
time is less than the specified time. The noted numbers of defects 
are plotted on the control chart to decide the state of the process. 
Previous studies [6] have proposed an attribute chart for Weibull 
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To detect indeterminacy effect in the manufacturing process, attribute control chart using neutrosophic Weibull distribution is 
proposed in this paper. To make the attribute control chart more efficient for persistent shifts in the industrial process, an attribute 
control chart using Weibull distribution has been proposed recently. In this study, a neutrosophic Weibull distribution-based 
attribute control chart develop for efficient monitoring of the process. The indeterminacy effect was studied with the control 
chart’s performance using characteristics of run length. In addition, the proposed chart effectively detected shifts in uncertainty. 
The relative efficiency of the proposed structure is compared with the existing attribute control chart under the Weibull time- 
truncated life test. The relative analysis reveals that the proposed time-truncated control chart for Weibull distribution under 
uncertainty design performance more efficiently than the existing counterparts. From the comparison, the proposed chart pro- 
vides smaller values for the out-of-control average run length as compared to the existing attribute control chart. An illustrative 
application related to automobile manufacturing is also incorporated to demonstrate the proposal. 


© 2021 The Authors. Published by Atlantis Press B.V. 


This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/). 


distribution. Specifically, some studies [7] designed control charts 
using censored data. Others [8] have studied time-truncated charts 
for exponentiated half logistic distribution and Dagum distribu- 
tion, respectively. Further, [9-12] show further control chart appli- 
cations for various statistical distributions. 


The Shewhart control charts given in the literature are unable to 
apply if some observations in the production data are uncertain or 
fuzzy. Fuzzy-based control charts are applied when uncertainty is 
presented in observations and parameters. According to [13], “fuzzy 
control charts are more sensitive than traditional ones; hence, they 
provide better quality products.” Furthermore, Darestani et al. [14] 
proposed a u-chart using a fuzzy approach, whereas [15] proposed 
a chart using fuzzy logic and [16] proposed a c-chart using fuzzy 
logic [17]. More work can be seen in [17-21]. 


The major drawback of a fuzzy-based control chart is that it is less 
informative than other charts. As such, Smarandache [22] men- 
tioned that the generalized form of fuzzy logic is known as neutro- 
sophic logic. Moreover, Smarandache and Khalid [23] stated that 
neutrosophic logic is more efficient than fuzzy logic and interval- 
based analysis. Other researchers [24] have discussed neutrosophic 
logic applications. Neutrosophic statistics is a branch of statistics 
that analyzes neutrosophic data, and many studies [25] have pre- 
sented methods to analyze the indeterminacy data. Other works 
[26] have introduced indeterminacy in an attribute chart and pro- 
posed attribute charts for neutrosophic statistics [27]. Moreover, 
Khan et al. [28] proposed an S-chart using neutrosophic statistics. 
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Lastly, Aslam et al. [29] proposed an efficient attribute chart. More 
information in neutrosophic statistics can be seen in [30,31]. 


The attribute control charts under neutrosophic statistics are 
presented in the literature. No previous work has designed a 
time-truncated attribute control chart for neutrosophic Weibull 
distribution using neutrosophic statistics. In this paper, we 
designed a time-truncated attribute control chart using neutro- 
sophic statistics. The effect of the indeterminacy parameter was 
studied on the average run length (ARL). The advantages of 
the proposed chart are also discussed. In next section, we describe 
the design of the proposed chart. In Section 3, the comparison of 
the proposed chart with the counterpart chart and simulation study 
is described. The Illustrative example is described in Section 4. In 
the final section, the conclusion and future recommendations are 
displayed. 


2. DESIGN OF THE PROPOSED CHART 


First, we introduce the neutrosophic Weibull distribution. Then, we 
present the design of the current chart. 


Suppose that x,e [xr xul be a neutrosophic random variable 
having neutrosophic Weibull distribution, where x; and xy are 
the lower and upper values of neutrosophic random variable. 
The Weibull distribution neutrosophic probability density func- 
tion (NPDF) and neutrosophic cumulative distribution function 
(NCDF) are defined as follows 


ren) = {(2) (2 EN} of) YE} 
Iye [rs Iy] sn [xr xy] 


(1) 


F (xy) =1- f (1 +1y)} 


+y; Ine |I Iy] »Xye [xr xy] 


(2) 


where p is the shape parameter, a is the scale parameter, and Iy 
shows the neutrosophic interval. The neutrosophic mean, say 4y of 
the neutrosophic Weibull distribution is given by 


Hy = a (1 + 1/8) (1+ Ty) (3) 


whereT (.) denotes the gamma function. Let 4yo target the product's 
lifetime and xy be a truncated time. The probability p y that an item 
fails by time xy is given as 


xy \F 
ies je) (1 +1) } ie ie tal asain we 
(4) 


Let Xyo = Gftyp be the termination time, where a is constant and we 
express the unknown scale parameter « in term of uy using Eq. (3). 
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Then, Eq. (4) can be rewritten as 


Pad ( fa.) ranyon) ( 


Py=l- 1+Iy) 


(5) 
+y; Ine [I Ty] » XNE [xr xy] 
When fy = Hyo then the probability in Eq. (5) reduces to 
2 {er (Eao 141 \ 
Pro ( N) (6) 
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The np attribute control chart having lower control limit (LCL) and 
upper control limit (UCL) for the neutrosophic Weibull distribution 
is explained in the following steps: 


Step 1. Count defective items of the products (D) from the selected 
sample. 


Step 2. If D > UCL or D < LCL, the process is considered out of 
control. The process is considered in control if LCL < D < UCL. 


Note that D follows a neutrosophic binomial distribution with 
parameters n and in-control probability py. Therefore, the control 
limits of the proposed np control chart are as follows: 


UCL = npyo + ky/npyo (1 -= Pro) (7a) 


LCL = max fo. NPno — KAS nPNo (1 = Px) (7b) 


where k is the coefficient or control constant of the control limits 
to be resolute. The pyy values are normally unknown; therefore, the 
averages of failed items (D = defective items) such as Dare taken. 
Thus, the neutrosophic control limits for practical use are 


UCLy =D + ky/D(1- D/n) (8a) 
LCLy = max 0,.D-kyD(1- 2) (8b) 


The probability for the in-control process (i.e., P.,,,) is given by 


P? „ =P (LCL < D < UCL|pyo) 
UCL n n-d (9) 
= Da ( 7) ol -= Pro) 


If LCL = 0, d should be 0 when applied. Let yy, be the shifted 
average and the probability at the new mean is evaluated by 


af ( He E b 
Pm =1— ye (me) eae (1+Iy) (10) 


+Iy3 Ine |I Iy] » XyE [xr xy] 
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If the shifted mean is 4y; = flmo for a constant f, then Eq. (10) can 
be rewritten as 

Pm =l- {PO Tao) (1 + Iy)} (11) 


+Iy; Ive [IIg] , XNE [xr xy] 


The in-process probability for the shifted mean can be written as 


Pi, =P (LCL < D < UCL|py;) 








Nin 
UCL n\ 4 fed (12) 
= Daa ( d ) PwC — Pm) 
The in-control ARL can be calculated as 
1 
ARL; = (13) 
1— Prin 
The ARL for the shifted process can be calculated as 
ARL, = —} (14) 
! l- Piin 


ARL values for various values (e.g., Iy, f, and p) are shown in 
Tables 1-3, where the following behavior can be noted. 


1. When the indeterminacy parameter Iy increases, the decreas- 
ing trend in ARL, is noted. 

2. For the same values, the values of ARL, decreases as p 
increases. 

The values of k, ARL}, and ARL, can be obtained through the fol- 

lowing steps: 


Step 1: Fix the values of p, n, and Iy, thus determining the values of 
k via the grid search method. 


Step 2: Several values of k are noted during simulation, where 
ARL, = 1 is a specified ARL value. 


Step 3: Choose the k value when ARL; is close to rọ. 


Step 4: Determine ARL, values for various values of f. 


Table 1 ARLs of the proposed chart when ry = 370. 
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3. COMPARATIVE STUDIES 


Our proposed control chart becomes the same chart proposed in 
[6] when Iy = 0. Table 4 presents both control charts when Iy = 
0.1 and Iy, = 0.5. Table 4 shows that the proposed control chart 
provides smaller ARLs values when compared to the chart pro- 
posed by [6]. For example, when f = 0.9, p = 1.1, n = 30, and 
Iy = 0.1, the proposed control chart detects the first out-of-control 
value at the 151st sample, whereas the existing chart proposed by 
[6] detects the first out-of-control value at the 177th sample. Simi- 
larly, when f = 0.9, p = 1.1, n = 30, and Iy = 0.5, the proposed 
control chart detects the first out-of-control value at the 131st sam- 
ple, whereas and the existing chart proposed by [6] detects the first 
out-of-control value at the 177th sample. From this study, it can be 
seen that the proposed chart provides a quick indication about the 
shift in the process when compared to the control chart proposed 
by [6]. This study thus proved the efficiency of the proposed chart. 














3.1. Simulation Study 


The performance of the proposed control chart was compared with 
the existing chart proposed by [6] using simulated data. The data 
was generated when a = 0.1548, n = 30,8 = 1, and Iy = 
0.1. Among 50 observations, the first 20 values were generated when 
the process was in-control (IC) state when fy = 1. Further, the 
next 30 observations were generated from the shifted process when 
f= 0.8 by using the R software. The values of the number of defec- 
tive items (D) under the proposed scheme are plotted in Figure 1 
and the existing control chart displayed in Figure 2. Figure 1 indi- 
cates the shift at the 31st sample; the existing chart proposed by [6] 
indicates no shift. By comparing both figures, we can conclude that 
the proposed control chart indicates an issue in the existing con- 
trol chart. Therefore, the application of the proposed control chart 
under uncertainty is helpful for minimizing the number of non- 
conforming items and got conforming products in less time. 


4. ILLUSTRATIVE EXAMPLE 


The proposed control chart was applied in an automobile 
manufacturing company located in South Korea. The company was 
interested in monitoring the service time in months for specific sub- 
systems [6]. The same service data was applied by [6], who showed 














B=1,n=30 
k = 3.26994, k = 3.12394, k = 3.23313, k = 3.22172, 
f Ty =0.1, a = 0.1545 Ty = 0.2, a = 0.1289 Ty = 0.4, a = 0.1082 Ty = 0.5, a = 0.09385 
ARL 

1.0 374.93 373.97 373.41 370.43 
0.9 165.38 163.88 154.50 152.48 
0.8 70.06 68.93 61.44 60.33 

0.7 28.65 27.97 23.68 23.13 

0.6 11.47 11.11 9.03 8.79 

0.5 4.65 4.48 3.59 3.48 

0.4 2.07 1.99 1.65 1.61 

0.3 1.18 1.15 1.06 1.05 

0.2 1.00 1.00 1.00 1.00 

0.1 1.00 1.00 1.00 1.00 
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Table 2 ARLs of the proposed chart when ry = 370. 
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B=1.1,n=30 
k = 3.13684, Iy k = 3.19848, Iy k = 3.1863, Iy k = 3.21308, Iy 
f = 0.1, a = 0.1883 = 0.2, a = 0.1584 = 0.4, a = 0.1332 = 0.5, a = 0.1310 
ARL 

1.0 371.62 371.68 371.37 374.92 

0.9 151.50 150.41 141.13 131.00 

0.8 59.51 58.61 51.77 46.61 

0.7 22.75 22.22 18.61 16.04 

0.6 8.67 8.40 6.78 5.68 

0.5 3.47 3.35 2.69 2.27 

0.4 1.64 1.50 1.34 1.21 

0.3 1.07 1.00 1.01 1.00 

0.2 1.00 1.00 1.00 1.00 

0.1 1.00 1.00 1.00 1.00 

Table3 ARLs of the proposed chart when rọ = 370. 
p=2,n=30 
k = 3.58595, Iy k = 3.58132 Iy k = 2.89913, Iy k = 3.3222, Iy 
f = 0.1, a = 0.2167 = 0.2, a = 0.2224 = 0.4, a = 0.2472 = 0.5, a = 0.2428 
ARL 

1.0 372.83 370.03 371.83 370.58 

0.9 136.70 120.64 88.27 79.58 

0.8 47.76 37.90 21.24 17.70 

0.7 16.20 11.85 5.64 4.53 

0.6 5.58 3.96 1.94 1.62 

0.5 2.17 1.64 1.09 1.03 

0.4 1.17 1.05 1.00 1.00 

0.3 1.00 1.00 1.00 1.001 

0.2 1.00 1.00 1.00 1.00 

0.1 1.00 1.00 1.00 1.00 

Table 4 Comparison of the proposed chart with the existing chart in ARLs. 
[6] chart ( Iy = 0) The proposed chart 
p=1 p=1.1 Iy=0.1 Iy=0.5 
f k = 3.40009, k = 3.33227, k = 3.26994, Iy k = 3.13684, f k = 3.22172, B k = 3.21308, f 
a= 0.1344 a= 0.1394 = 0.1, a= 0.1545 =1.1, a= 0.1883 = 1, a=0.09385 =1.1,a=0.1310 
ARL 

1.0 370.32 370.80 374.93 371.62 370.43 374.92 
0.9 180.12 177.71 165.38 151.50 152.48 131.00 
0.8 83.88 81.21 70.06 59.51 60.33 46.61 
0.7 37.40 35.50 28.65 22.75 23.13 16.04 
0.6 16.05 14.94 11.47 8.67 8.79 5.68 
0.5 6.74 6.18 4.65 3.47 3.48 2.27 
0.4 2.91 2.66 2.07 1.64 1.61 1.21 
0.3 1.45 1.35 1.18 1.07 1.05 1.00 
0.2 1.03 1.01 1.00 1.00 1.00 1.00 
0.1 1.00 1.00 1.00 1.00 1.00 1.00 





that the data followed the Weibull distribution with J = 1. For this 
study, let a = 0.1883, k = 3.1368, n = 30, py = 0.1722, and 
Iy = 0.1, which leads to the value of 6 = 1(14+0.1) = 1.1 for 
the neutrosophic Weibull distribution. Let the value of zy) = 60, 
which leads to time Xpo = Gfyy = 0.1883 * 60 = 11.29. From 400 


observations, the statistic D is made by counting the observation as 
a failure if it is smaller than 11.29 months. The numbers of defec- 
tives D are shown in Figures 3 and 4, and also listed as follows: 


3, 1, 4, 4, 3, 3, 3, 6, 4, 3, 5, 2, 3, 3, 4, 4, 4, 4, 6, 1, 4, 5, 5, 6, 7, 7, 2, 1,7, 
8, 3, 3, 3, 9, 4, 5, 10, 5, 5, 7. 
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Figure 1 


The proposed chart for simulation data. 
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Figure 2 The [6] chart for simulation data. 


The control limits for the proposed chart are UCL = 8 and LCL = 
0 when the values of parameters values are used from Table 2. 
The proposed chart is displayed in Figure 3 and the existing chart 
is shown in Figures 4, respectively. Figure 3 shows the shift in 
service time means the shift values detect out of control at 31st 
sample. The existing chart indicates that no action is needed for 
service time because all values are in control. By comparing both 
charts, it can be seen that the proposed control chart indicates an 
issue in monthly service time for the specific subsystem of the car 
while the existing chart does not show any issue in monthly service 
time for the specific subsystem of the car. 


5. CONCLUSIONS 


The time-truncated attribute control chart was herein presented 
using the neutrosophic Weibull distribution. The measurement’s 


indeterminacy effect was studied on the performance of the pro- 
posed control chart with the counterpart control chart. The pro- 
posed chart was the generalized version of the attribute control 
chart under a neutrosophic environment. From the results pre- 
sented in Tables 1-4, we observed that the indeterminacy parame- 
ter significantly affected ARL values. The out-of-control ARL val- 
ues reduced when indeterminacy increased. The comparative study 
proved the current chart’s efficiency. The proposed control chart 
was proven to monitor service time in the automobile industry. The 
proposed control chart has limitations, i.e., it can only be applied 
when the life/service time follows the neutrosophic Weibull dis- 
tribution. Secondly, the proposed control chart cannot be applied 
for variable data. The proposed control chart can be applied in 
the automobile industry, aircraft industry, and mobile industry 
for monitoring the defective items. In the future, the proposed 
control chart should be used for other neutrosophic statistical dis- 
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Figure 3 The proposed chart for real data. 
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Figure 4 The existing chart for real data. 


tributions. Following [7], the proposed control chart used an expo- 
nentially weighted moving average (EWMA) statistic and cumula- 
tive sum (CUSUM) statistic, both of which should be considered in 
future research. 
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